\(\int \frac {\sinh ^3(a+b x^2)}{x^3} \, dx\) [21]

Optimal result

Integrand size = 14, antiderivative size = 91 \[ \int \frac {\sinh ^3\left (a+b x^2\right )}{x^3} \, dx=-\frac {3}{8} b \cosh (a) \text {Chi}\left (b x^2\right )+\frac {3}{8} b \cosh (3 a) \text {Chi}\left (3 b x^2\right )+\frac {3 \sinh \left (a+b x^2\right )}{8 x^2}-\frac {\sinh \left (3 \left (a+b x^2\right )\right )}{8 x^2}-\frac {3}{8} b \sinh (a) \text {Shi}\left (b x^2\right )+\frac {3}{8} b \sinh (3 a) \text {Shi}\left (3 b x^2\right ) \]

[Out]

-3/8*b*Chi(b*x^2)*cosh(a)+3/8*b*Chi(3*b*x^2)*cosh(3*a)-3/8*b*Shi(b*x^2)*sinh(a)+3/8*b*Shi(3*b*x^2)*sinh(3*a)+3
/8*sinh(b*x^2+a)/x^2-1/8*sinh(3*b*x^2+3*a)/x^2

Rubi [A] (verified)

Time = 0.14 (sec) , antiderivative size = 91, normalized size of antiderivative = 1.00, number of steps used = 12, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.429, Rules used = {5448, 5428, 3378, 3384, 3379, 3382} \[ \int \frac {\sinh ^3\left (a+b x^2\right )}{x^3} \, dx=-\frac {3}{8} b \cosh (a) \text {Chi}\left (b x^2\right )+\frac {3}{8} b \cosh (3 a) \text {Chi}\left (3 b x^2\right )-\frac {3}{8} b \sinh (a) \text {Shi}\left (b x^2\right )+\frac {3}{8} b \sinh (3 a) \text {Shi}\left (3 b x^2\right )+\frac {3 \sinh \left (a+b x^2\right )}{8 x^2}-\frac {\sinh \left (3 \left (a+b x^2\right )\right )}{8 x^2} \]

[In]

Int[Sinh[a + b*x^2]^3/x^3,x]

[Out]

(-3*b*Cosh[a]*CoshIntegral[b*x^2])/8 + (3*b*Cosh[3*a]*CoshIntegral[3*b*x^2])/8 + (3*Sinh[a + b*x^2])/(8*x^2) -
Sinh[3*(a + b*x^2)]/(8*x^2) - (3*b*Sinh[a]*SinhIntegral[b*x^2])/8 + (3*b*Sinh[3*a]*SinhIntegral[3*b*x^2])/8

Rule 3378

Int[((c_.) + (d_.)*(x_))^(m_)*sin[(e_.) + (f_.)*(x_)], x_Symbol] :> Simp[(c + d*x)^(m + 1)*(Sin[e + f*x]/(d*(m
+ 1))), x] - Dist[f/(d*(m + 1)), Int[(c + d*x)^(m + 1)*Cos[e + f*x], x], x] /; FreeQ[{c, d, e, f}, x] && LtQ[
m, -1]

Rule 3379

Int[sin[(e_.) + (Complex[0, fz_])*(f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[I*(SinhIntegral[c*f*(fz/
d) + f*fz*x]/d), x] /; FreeQ[{c, d, e, f, fz}, x] && EqQ[d*e - c*f*fz*I, 0]

Rule 3382

Int[sin[(e_.) + (Complex[0, fz_])*(f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[CoshIntegral[c*f*(fz/d)
+ f*fz*x]/d, x] /; FreeQ[{c, d, e, f, fz}, x] && EqQ[d*(e - Pi/2) - c*f*fz*I, 0]

Rule 3384

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Dist[Cos[(d*e - c*f)/d], Int[Sin[c*(f/d) + f*x]
/(c + d*x), x], x] + Dist[Sin[(d*e - c*f)/d], Int[Cos[c*(f/d) + f*x]/(c + d*x), x], x] /; FreeQ[{c, d, e, f},
x] && NeQ[d*e - c*f, 0]

Rule 5428

Int[(x_)^(m_.)*((a_.) + (b_.)*Sinh[(c_.) + (d_.)*(x_)^(n_)])^(p_.), x_Symbol] :> Dist[1/n, Subst[Int[x^(Simpli
fy[(m + 1)/n] - 1)*(a + b*Sinh[c + d*x])^p, x], x, x^n], x] /; FreeQ[{a, b, c, d, m, n, p}, x] && IntegerQ[Sim
plify[(m + 1)/n]] && (EqQ[p, 1] || EqQ[m, n - 1] || (IntegerQ[p] && GtQ[Simplify[(m + 1)/n], 0]))

Rule 5448

Int[((e_.)*(x_))^(m_.)*((a_.) + (b_.)*Sinh[(c_.) + (d_.)*(x_)^(n_)])^(p_), x_Symbol] :> Int[ExpandTrigReduce[(
e*x)^m, (a + b*Sinh[c + d*x^n])^p, x], x] /; FreeQ[{a, b, c, d, e, m}, x] && IGtQ[p, 1] && IGtQ[n, 0]

Rubi steps \begin{align*} \text {integral}& = \int \left (-\frac {3 \sinh \left (a+b x^2\right )}{4 x^3}+\frac {\sinh \left (3 a+3 b x^2\right )}{4 x^3}\right ) \, dx \\ & = \frac {1}{4} \int \frac {\sinh \left (3 a+3 b x^2\right )}{x^3} \, dx-\frac {3}{4} \int \frac {\sinh \left (a+b x^2\right )}{x^3} \, dx \\ & = \frac {1}{8} \text {Subst}\left (\int \frac {\sinh (3 a+3 b x)}{x^2} \, dx,x,x^2\right )-\frac {3}{8} \text {Subst}\left (\int \frac {\sinh (a+b x)}{x^2} \, dx,x,x^2\right ) \\ & = \frac {3 \sinh \left (a+b x^2\right )}{8 x^2}-\frac {\sinh \left (3 \left (a+b x^2\right )\right )}{8 x^2}-\frac {1}{8} (3 b) \text {Subst}\left (\int \frac {\cosh (a+b x)}{x} \, dx,x,x^2\right )+\frac {1}{8} (3 b) \text {Subst}\left (\int \frac {\cosh (3 a+3 b x)}{x} \, dx,x,x^2\right ) \\ & = \frac {3 \sinh \left (a+b x^2\right )}{8 x^2}-\frac {\sinh \left (3 \left (a+b x^2\right )\right )}{8 x^2}-\frac {1}{8} (3 b \cosh (a)) \text {Subst}\left (\int \frac {\cosh (b x)}{x} \, dx,x,x^2\right )+\frac {1}{8} (3 b \cosh (3 a)) \text {Subst}\left (\int \frac {\cosh (3 b x)}{x} \, dx,x,x^2\right )-\frac {1}{8} (3 b \sinh (a)) \text {Subst}\left (\int \frac {\sinh (b x)}{x} \, dx,x,x^2\right )+\frac {1}{8} (3 b \sinh (3 a)) \text {Subst}\left (\int \frac {\sinh (3 b x)}{x} \, dx,x,x^2\right ) \\ & = -\frac {3}{8} b \cosh (a) \text {Chi}\left (b x^2\right )+\frac {3}{8} b \cosh (3 a) \text {Chi}\left (3 b x^2\right )+\frac {3 \sinh \left (a+b x^2\right )}{8 x^2}-\frac {\sinh \left (3 \left (a+b x^2\right )\right )}{8 x^2}-\frac {3}{8} b \sinh (a) \text {Shi}\left (b x^2\right )+\frac {3}{8} b \sinh (3 a) \text {Shi}\left (3 b x^2\right ) \\ \end{align*}

Mathematica [A] (verified)

Time = 0.08 (sec) , antiderivative size = 90, normalized size of antiderivative = 0.99 \[ \int \frac {\sinh ^3\left (a+b x^2\right )}{x^3} \, dx=-\frac {3 b x^2 \cosh (a) \text {Chi}\left (b x^2\right )-3 b x^2 \cosh (3 a) \text {Chi}\left (3 b x^2\right )-3 \sinh \left (a+b x^2\right )+\sinh \left (3 \left (a+b x^2\right )\right )+3 b x^2 \sinh (a) \text {Shi}\left (b x^2\right )-3 b x^2 \sinh (3 a) \text {Shi}\left (3 b x^2\right )}{8 x^2} \]

[In]

Integrate[Sinh[a + b*x^2]^3/x^3,x]

[Out]

-1/8*(3*b*x^2*Cosh[a]*CoshIntegral[b*x^2] - 3*b*x^2*Cosh[3*a]*CoshIntegral[3*b*x^2] - 3*Sinh[a + b*x^2] + Sinh
[3*(a + b*x^2)] + 3*b*x^2*Sinh[a]*SinhIntegral[b*x^2] - 3*b*x^2*Sinh[3*a]*SinhIntegral[3*b*x^2])/x^2

Maple [A] (verified)

Time = 0.62 (sec) , antiderivative size = 123, normalized size of antiderivative = 1.35

method	result	size
risch	\(\frac {-3 \,{\mathrm e}^{-3 a} \operatorname {Ei}_{1}\left (3 x^{2} b \right ) b \,x^{2}+3 \,\operatorname {Ei}_{1}\left (x^{2} b \right ) {\mathrm e}^{-a} b \,x^{2}+3 \,{\mathrm e}^{a} \operatorname {Ei}_{1}\left (-x^{2} b \right ) b \,x^{2}-3 \,{\mathrm e}^{3 a} \operatorname {Ei}_{1}\left (-3 x^{2} b \right ) b \,x^{2}+{\mathrm e}^{-3 x^{2} b -3 a}-3 \,{\mathrm e}^{-x^{2} b -a}+3 \,{\mathrm e}^{x^{2} b +a}-{\mathrm e}^{3 x^{2} b +3 a}}{16 x^{2}}\)	\(123\)

[In]

int(sinh(b*x^2+a)^3/x^3,x,method=_RETURNVERBOSE)

[Out]

1/16*(-3*exp(-3*a)*Ei(1,3*x^2*b)*b*x^2+3*Ei(1,x^2*b)*exp(-a)*b*x^2+3*exp(a)*Ei(1,-x^2*b)*b*x^2-3*exp(3*a)*Ei(1
,-3*x^2*b)*b*x^2+exp(-3*b*x^2-3*a)-3*exp(-b*x^2-a)+3*exp(b*x^2+a)-exp(3*b*x^2+3*a))/x^2

Fricas [A] (verification not implemented)

none

Time = 0.25 (sec) , antiderivative size = 160, normalized size of antiderivative = 1.76 \[ \int \frac {\sinh ^3\left (a+b x^2\right )}{x^3} \, dx=-\frac {2 \, \sinh \left (b x^{2} + a\right )^{3} - 3 \, {\left (b x^{2} {\rm Ei}\left (3 \, b x^{2}\right ) + b x^{2} {\rm Ei}\left (-3 \, b x^{2}\right )\right )} \cosh \left (3 \, a\right ) + 3 \, {\left (b x^{2} {\rm Ei}\left (b x^{2}\right ) + b x^{2} {\rm Ei}\left (-b x^{2}\right )\right )} \cosh \left (a\right ) + 6 \, {\left (\cosh \left (b x^{2} + a\right )^{2} - 1\right )} \sinh \left (b x^{2} + a\right ) - 3 \, {\left (b x^{2} {\rm Ei}\left (3 \, b x^{2}\right ) - b x^{2} {\rm Ei}\left (-3 \, b x^{2}\right )\right )} \sinh \left (3 \, a\right ) + 3 \, {\left (b x^{2} {\rm Ei}\left (b x^{2}\right ) - b x^{2} {\rm Ei}\left (-b x^{2}\right )\right )} \sinh \left (a\right )}{16 \, x^{2}} \]

[In]

integrate(sinh(b*x^2+a)^3/x^3,x, algorithm="fricas")

[Out]

-1/16*(2*sinh(b*x^2 + a)^3 - 3*(b*x^2*Ei(3*b*x^2) + b*x^2*Ei(-3*b*x^2))*cosh(3*a) + 3*(b*x^2*Ei(b*x^2) + b*x^2
*Ei(-b*x^2))*cosh(a) + 6*(cosh(b*x^2 + a)^2 - 1)*sinh(b*x^2 + a) - 3*(b*x^2*Ei(3*b*x^2) - b*x^2*Ei(-3*b*x^2))*
sinh(3*a) + 3*(b*x^2*Ei(b*x^2) - b*x^2*Ei(-b*x^2))*sinh(a))/x^2

Sympy [F]

\[ \int \frac {\sinh ^3\left (a+b x^2\right )}{x^3} \, dx=\int \frac {\sinh ^{3}{\left (a + b x^{2} \right )}}{x^{3}}\, dx \]

[In]

integrate(sinh(b*x**2+a)**3/x**3,x)

[Out]

Integral(sinh(a + b*x**2)**3/x**3, x)

Maxima [A] (verification not implemented)

none

Time = 0.29 (sec) , antiderivative size = 58, normalized size of antiderivative = 0.64 \[ \int \frac {\sinh ^3\left (a+b x^2\right )}{x^3} \, dx=\frac {3}{16} \, b e^{\left (-3 \, a\right )} \Gamma \left (-1, 3 \, b x^{2}\right ) - \frac {3}{16} \, b e^{\left (-a\right )} \Gamma \left (-1, b x^{2}\right ) - \frac {3}{16} \, b e^{a} \Gamma \left (-1, -b x^{2}\right ) + \frac {3}{16} \, b e^{\left (3 \, a\right )} \Gamma \left (-1, -3 \, b x^{2}\right ) \]

[In]

integrate(sinh(b*x^2+a)^3/x^3,x, algorithm="maxima")

[Out]

3/16*b*e^(-3*a)*gamma(-1, 3*b*x^2) - 3/16*b*e^(-a)*gamma(-1, b*x^2) - 3/16*b*e^a*gamma(-1, -b*x^2) + 3/16*b*e^
(3*a)*gamma(-1, -3*b*x^2)

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 223 vs. \(2 (80) = 160\).

Time = 0.28 (sec) , antiderivative size = 223, normalized size of antiderivative = 2.45 \[ \int \frac {\sinh ^3\left (a+b x^2\right )}{x^3} \, dx=\frac {3 \, {\left (b x^{2} + a\right )} b^{2} {\rm Ei}\left (3 \, b x^{2}\right ) e^{\left (3 \, a\right )} - 3 \, a b^{2} {\rm Ei}\left (3 \, b x^{2}\right ) e^{\left (3 \, a\right )} - 3 \, {\left (b x^{2} + a\right )} b^{2} {\rm Ei}\left (-b x^{2}\right ) e^{\left (-a\right )} + 3 \, a b^{2} {\rm Ei}\left (-b x^{2}\right ) e^{\left (-a\right )} + 3 \, {\left (b x^{2} + a\right )} b^{2} {\rm Ei}\left (-3 \, b x^{2}\right ) e^{\left (-3 \, a\right )} - 3 \, a b^{2} {\rm Ei}\left (-3 \, b x^{2}\right ) e^{\left (-3 \, a\right )} - 3 \, {\left (b x^{2} + a\right )} b^{2} {\rm Ei}\left (b x^{2}\right ) e^{a} + 3 \, a b^{2} {\rm Ei}\left (b x^{2}\right ) e^{a} - b^{2} e^{\left (3 \, b x^{2} + 3 \, a\right )} + 3 \, b^{2} e^{\left (b x^{2} + a\right )} - 3 \, b^{2} e^{\left (-b x^{2} - a\right )} + b^{2} e^{\left (-3 \, b x^{2} - 3 \, a\right )}}{16 \, b^{2} x^{2}} \]

[In]

integrate(sinh(b*x^2+a)^3/x^3,x, algorithm="giac")

[Out]

1/16*(3*(b*x^2 + a)*b^2*Ei(3*b*x^2)*e^(3*a) - 3*a*b^2*Ei(3*b*x^2)*e^(3*a) - 3*(b*x^2 + a)*b^2*Ei(-b*x^2)*e^(-a
) + 3*a*b^2*Ei(-b*x^2)*e^(-a) + 3*(b*x^2 + a)*b^2*Ei(-3*b*x^2)*e^(-3*a) - 3*a*b^2*Ei(-3*b*x^2)*e^(-3*a) - 3*(b
*x^2 + a)*b^2*Ei(b*x^2)*e^a + 3*a*b^2*Ei(b*x^2)*e^a - b^2*e^(3*b*x^2 + 3*a) + 3*b^2*e^(b*x^2 + a) - 3*b^2*e^(-
b*x^2 - a) + b^2*e^(-3*b*x^2 - 3*a))/(b^2*x^2)

Mupad [F(-1)]

Timed out. \[ \int \frac {\sinh ^3\left (a+b x^2\right )}{x^3} \, dx=\int \frac {{\mathrm {sinh}\left (b\,x^2+a\right )}^3}{x^3} \,d x \]

[In]

int(sinh(a + b*x^2)^3/x^3,x)

[Out]

int(sinh(a + b*x^2)^3/x^3, x)